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Abstract. 

High-temperature cuprate superconductors have been known to exhibit significant pressure 
effects. In order to fathom the origin of why and how T c is affected by pressure, we have recently 
studied the pressure effects on T c adopting a model that contains two copper d-orbitals derived 
from first-principles band calculations, where the d z i orbital is considered on top of the usually 
considered d x i_ y 2 orbital. In that paper, we have identified two origins for the T c enhancement 
under hydrostatic pressure: (i) while at ambient pressure the smaller the hybridization of other 
orbital components the higher the T c , an application of pressure acts to reduce the multiorbital 
mixing on the Fermi surface, which we call the orbital distillation effect, and (ii) the increase 
of the band width with pressure also contributes to the enhancement. In the present paper, 
we further elaborate the two points. As for point (i), while the reduction of the apical oxygen 
height under pressure tends to increase the d z 2 mixture, hence to lower T c , here we show that 
this effect is strongly reduced in bi-layer materials due to the pyramidal coordination of oxygen 
atoms. As for point (ii), we show that the enhancement of T c due to the increase in the band 
width is caused by the effect that the many-body renormalization arising from the self-energy 
is reduced. 



1. Introduction 

Although many kinds of superconductors have been discovered, the superconducting transition 
temperature T c of the cuprate superconductors still remains to be the highest, and the possibility 
of further enhancing T c still attracts much attention. To achieve higher T c , it is important to 
understand the key parameters that control T c , and from this viewpoint, there have been studies 
that have focused on the correlation between T c and the lattice parameters such as the in-plain 
bond length(7)[H E] or the apical oxygen height(/i ) [3], SI El d El [10l [11] measured from 
the Cu02 planes. Theoretically, we have studied the material dependence of T c in refs.fTD] ITT], 
and introduced a two-orbital model that takes into account the d x 2_ y 2 and d z 2 Wannier orbitals. 
In most of the theories of the cuprates, only the d x 2_ y 2 (and the hybridized oxygen p) orbital 



is considered, but actually it has been noticed from the early days that in (La,Sr)2Cu04, which 
has small ho and relatively low T c , there is a strong mixture of the d z i orbital component near 
the Fermi level[12j [T3"l 114] . In fact, nowadays there have been more studies that focus on the 
d z 2 orbital or the related apical oxygen El H51 [THl El US]- In refs.JTOl E], we have shown 
that this mixture of the d z i orbital component near the Fermi level works destructively against 
d— wave superconductivity, and hence this is the main reason of the material dependence of T c . 

Studying the pressure effect on T c is an in situ way of attacking the problem of the correlation 
between T c and the lattice structure. It is well known that in most of the cuprates, application 
of pressure enhances T^ |19l I2U] . On the other hand, it has been revealed that the pressure effect 
exhibits strong anisotropy[2TJ [221 I n our recent theoretical study on the pressure effect 
in the cuprates, it has been revealed that in addition to the d z i effect, the roundness of the 
Fermi surface, which is controlled by Cu-As-d x 2_ y 2 hybridization, and also the band width are 
important parameters that governs T c under pressure[24j. 

In the present paper, after briefly reviewing the main results obtained in ref , |24j . we 
study more closely the effect of the band width by applying the fluctuation exchange 
approximation(FLEX) [25"1 [26] [27] to the two-orbital model. Secondly, we will explain the 
difference of the d z 2 orbital effect between the multi-layer and single-layer systems under 
hydrostatic pressure. 

2. Calculation Method 

2.1. Determination of the Crystal Structure Under Pressure 

We first obtain the crystal structure of the single-layer cuprates La2Cu04 and Hg2BaCu04 under 
ambient pressure. Namely, we calculate the total energy by first principles calculation [28 j varying 
the lattice constants, and fit the result by the standard Burch-Marnaghan formula[29j to obtain 
the structure at the most stable point. From such calculation, we can obtain the crystal structure 
within one percent discrepancy from the lattice constants determined experimentally [3(3 13 Xj . 
To obtain the lattice structure under hydrostatic pressure, we optimize the Poisson's ratio and 
atomic position by first principles reducing the cell volume to 95% or 90% of the lattice structure 
under ambient pressure. Since the compressibility in the cuprates is known to be about ~ 0.01 
GPa -1 [52]. V = 0.9Vo corresponds to a pressure of about 10 GPa. 

2.2. Construction of the two-orbital model and FLEX approximation 

Using the obtained crystal structure under ambient or hydrostatic pressure, we construct 
maximally localized Wannier orbitals[33"l [33] to extract the hopping parameters of the d x 2_ y 2- 
d z 2 two-orbital model[10|. As for the electron-electron interactions, we consider the on-site 
intraorbital Coulomb repulsion U, interorbital repulsion U', the Hund's coupling J and the 
pair-hopping J' . Here we also keep the orbital SU(2) requirement, U — U' = 2 J. We set U = 3.0 
eV, U' = 2.4 eV and J = J' = 0.3 eV unless mentioned otherwise. Estimates of U for the 
cuprates is 7 — lOt, t ~ 0.45eV (namely, U is about 3 ~ 4.5 eV), and J(J') — 0.1U, so the values 
adopted here are within the widely accepted range. 

We apply FLEX to this model to obtain the Green's function renormalized by the many-body 
self-energy correction. In FLEX, we define the spin and charge susceptibilities as follows; 
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here, are orbital indices. Considering the self-energy correction originating from these 

susceptibilities, we solve the Dyson equation in a self-consistent manner. Then the renormalized 
Green's function is substituted to the linearized Eliashberg equation for superconductivity, 



The maximum eigenvalue A in the above equation reaches unity at the superconducting 
transition temperature T = T c , so that A calculated at a fixed temperature can be used as a 
measure for T c . Here, we calculate A at T = 0.01 eV for La and T = 0.03 eV for Hg. The reason 
for this is because T c ~ 100 K of Hg is about three times larger than for La with T c ~ 40 K[35j. 

3. Results and Discussions: effect of band width 

3.1. Orbital Distillation 

In Figdl A is depicted as a function of the unit cell volume. A increases in La and in Hg 
cuprates with hydrostatic pressure, and this agrees with the well-known experimental result 
that T c goes up monotonically under pressure in the cuprates [TUJ ED]. To understand this result, 
we have introduced three important parameters, the level offset AE, the roundness of Fermi 
surface (r x 2_ y 2 defined below) and the band width T^JM]. To decompose the pressure effect 
into the contribution from the variation of these parameters, we vary each parameter "by hand" 
from the original value at V = Vq to the value at V = 0.9Vq separately, and calculate the 
variation in A. This result is also depicted in FigQ] as the length of arrows. AE is the on-site 
energy difference between d x 2_ y 2 and d z 2 Wannier orbitals in the two-orbital model, so this is 
a measure of the d z 2 orbital effect, which has been found to work destructively against d-wave 
superconductivity in our previous study [TO]. The effect of AE is dominant in La compound 
because La originally has small AE thereby suppressing T c , so its increase under pressure is 
effective for the T c enhancement. Note that in La system, AE increases under pressure because 
the absolute magnitude of the crystal field increases in spite of the reduction of ho /I. On 
the other hand, in the Hg compound the pressure effect through AE is negligible because AE 
is intrinsically large. Instead, the other two parameters are effective. r x 2_ y 2 is defined as 
r x 2_ y 2 = (|*2 1 + |*3|)/|*i|i where U is the i-th neighbor hopping within the d x 2_ y 2 orbital. The 
roundness of Fermi surface is enhanced by the increase of this value, and ti and £3 are mostly 
mediated by the Cu-4s orbital (which is effectively included in the two Wannier orbitals in the 
present model) with the path of d x 2_ y 2 — > 4s — > d x 2_ y 2[4\. It is known that the roundness 
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of the Fermi surface works against the spin-fluctuation-mediated superconductivity [36j, so the 
reduction of the 4s effect enhances T c . The 4s orbital effect is reduced by hydrostatic pressure 
because the energy level offset between the d x 2_ y 2 and the 4s orbital is enhanced when the 
oxygen ligands approach Cu. The effect of d z 2 and 4s orbitals put together, we can say that 
T c increases when the main band has more pure d x 2_ y 2 component, namely, when the "orbital 
distillation" takes place. 
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Figure 1. Eigenvalue A of the Eliashberg equation plotted against the volume compression 
V/Vq. Arrows indicate the contribution to the A variance from the parameters AE, r x i_ y 2 and 
W (see text). 



3.2. The Effect of Band Width W 

Now, we turn to the first main topic of this paper, i.e., the effect of the band width W. The 
band width W is defined as the energy range between the top and the bottom of the main 
band. It is evident that the band width is controlled by the in-plain bond length I, so that 
hydrostatic pressure enhances W. In Fig{TJ we can see that the increase in W results in an 
enhancement of A in Hg, while the opposite occurs for La. The reason for this can be understood 
as follows. Let us first start with the U dependence of T c for a fixed band width. In the top 
of FigJ^l we show the absolute value of the renormalized Green's function squared \G\ 2 in the 
Hg compound at (k,iu) = (ir, 0, mksT) for several values of 2 < U < 5 eV. The value of 
\G\ 2 monotonically decreases with larger U because the self-energy, which increases with U, 
suppresses \G\ 2 . On the other hand, the pairing interaction shown in the middle increases with 
U because the spin fluctuations develop monotonically. Consequently, y|G| 2 shown on the right 
exhibits a maximum at a certain U. Since F|G| 2 can be considered as a rough measure of the 
eigenvalue of the Eliashberg equation for d— wave superconductivity, A (and thus T c ) is expected 
to be maximized around a certain C/ max . 

In the left panels of FigJSl we show the U dependence of A for 3 < U < 5 eV for the two 
compounds with V = Vq and 0.9Vb- We can see that U < 5 eV lies on the left side of U max 
for La, while 3 eV< U is on the right side of U max in Hg. This means that larger values of 
U is necessary for La to be in the "strongly correlated regime" . This is because electrons can 
avoid the strong intraorbital repulsion within d x 2_ y 2 orbitals by using the d z 2 orbital degrees of 
freedom in materials with small AE. 

In the right panel of Figj3l we show the U dependence of the increment AA of the eigenvalue 
A induced by the increase of W under pressure. In this calculation, W is increased "by hand" 
up to its value at V = 0.9Vo, while the other two parameters are fixed at their original values. 
From this figure, we can see that the increase in W always enhances T c in the Hg cuprate within 
the realistic U range. For the La cuprate on the other hand, AA is negative for small values of 
U, and this is the reason why W affects superconductivity in opposite ways between La and Hg 
in FigdJ For larger values of U, however, AA becomes positive even for La. 
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Figure 2. U dependence of the FLEX calculation result for the Hg compound, (left) the 
absolute value of the renormalized Green's function squared \G\ 2 at (k,iu) = (ir,0,iirkBT), 
(middle) the effective pairing interaction V at (q,iu)) = (it,tt,0) and (right) the product 1/|G| 2 . 
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Figure 3. (Left) A against U for La (top) and Hg (bottom). Filled (opened) symbols displays 
the result for V = Vq(V = 0.9Vb). (Right) AA (increment in A when W is increased up to its 
value at V = 0.9Vq) against U. 



To see the effect of the band width more clearly, we provide in FigJ5]a schematic band width 
dependence of the U vs T c plot. Namely, T c should depend on U and W essentially in the form 
W f(U/W), where / is a certain function that gives the overall dependence of T c against the 
electron correlation strength. Therefore, as the pressure is applied, W increases so that C/ max 
(peak position of the curve) also increases accordingly keeping U ma _ x /W constant. At the same 
time, the absolute value of the maximized T c is enhanced by the application of the pressure 
because the entire energy scale increases (i.e., both U max and W are enhanced). 

4. The d z 2 Orbital Effect: Multi-layer vs Single-layer Systems 

So far, we have concentrated on single-layer systems. In this section, we will discuss the difference 
of the pressure effect through AE between single and multi-layer systems. Here, we focus on 
the comparison between single- and bi-layer La, Hg, Tl, Bi, and Y cuprates. We consider AEj, 
the energy difference between the d x 2_ y 2 and the d z 2 orbital in the d-p model which considers all 




Figure 4. Schematic figure of T c variation induced by applied pressure. Each line is a U vs. 
T c plot for fixed W, and W is larger for solid green — > dashed blue — > dotted red. 



of the Cu 3d and 2p orbitals explicitly. In the d-p model, the basis functions for the hopping 
part are close to the atomic orbitals, so AEd can be considered as the energy difference between 
the d x 2_ y 2 and the d z 2 atomic orbitals. AE is the energy difference between the d x 2_ y 2 and the 
d z 2 Wannier orbitals which effectively take the 02p orbitals into account, so AEd and AE are 
positively correlated 

In FigEl we show the relationship between the apical oxygen height ho and AEd- We can 
see that, while AEd is positively correlated with the apical oxygen height in both single- and 
bi-layer systems as expected, AEd is overall significantly greater in the bi-layer systems than in 
the single-layer systems. This is because bi-layer cuprates take pyramidal coordination of the 
oxygen ligands, while the single-layer cuprates take an octahedral one. Hence the effect of the 
apical oxygen should be weaker in the bi-layer systems, so that the "effective ho" is larger. As 
we have seen in the previous section, AE (and hence the d z 2 orbital) plays minor role in the 
T c enhancement under pressure in materials with large AE, so the effect of AE should become 
less relevant as the number of the layers increases. This can be considered as one reason why 
T c increases in spite of somewhat larger reduction in ho in bi-layer system than in single-layer 
system under hydrostatic pressure[T]. 




Figure 5. Relationship between the apical oxygen height /io[^J an d the level difference 
AEd [eV]. The squares(circles) show the value of AEd in each bi-layer(single-layer) system(the 
materials are indicated by arrows). The dashed(solid) line presents the calculation result when 
ho is varied hypothetically in bi-layer (single-layer) La system by hand. 



5. Conclusion 

In conclusion, we have performed first-principle calculation+FLEX study to analyze the pressure 
effect on T c of the cuprates. To explain the T c enhancement induced by hydrostatic pressure, we 
introduce three important parameters, level offset AE, roundness of the Fermi surface r x 2_ y 2 
and the band width W. AE and r x 2_ y 2 are the measure of "orbital purity" in the main band, 
and T c is enhanced by "orbital distillation" . We have shown that hydrostatic pressure enhances 
the distillation, thereby enhancing T c . We have also analyzed the effect of the band width W in 
detail. It has been shown that A (and hence T c ) is maximized around a certain U, which depends 
on the band width as well as AE. For materials with large AE such as the Hg compound, the 
increase of W under pressure leads to an enhancement of A for realistic values of U. 

We also discuss the difference between bi- or single-layer systems from the view point of AE. 
Even when ho is suppressed under pressure in bi-layer systems, the effect on T c should be small 
because of the pyramidal coordination of the oxygen ligands. 
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